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Abstract 

We study the the survival probability P{t) upto time t, of a test particle 
moving in a fluctuating external field. The particle moves according to 
some prescribed deterministic or stochastic rules and survives as long as 
the external field that it sees at its own location does not change sign. This 
is a natural generalization of the "static persistence" (when the particle is at 
rest) that has generated considerable recent interests. Two types of motions 
of the particle are considered. In one case, the particle adopts a strategy 
to live longer and in the other it just diffuses randomly. Three different 
external fields were considered: (i) the solution of diffusion equation, (ii) 
the "colour" profile of the g-state Potts model undergoing zero temperature 
coarsening dynamics and (iii) spatially uncorrelated Brownian signals. In 
most cases studied, P{t) ~ t^^" for large t. The exponent dm is calculated 
via numerically, analytically by approximate methods and in some cases 
exactly. It is shown in some special cases that the survival probability of 
the mobile particle is related to the persistence of special "patterns" present 
in the initial configuration of a phase ordering system. 
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I. Introduction 



Considerable interests have been generated recently in understanding the statistics 
of first passage events in spatially extended nonequilibrium systems. These sys- 
tems include the Ising or Potts model undergoing zero temperature phase order- 
ing dynamics |l|, |], simple diffusion equation with random initial conditions 0, 
H, several reaction-diffusion systems]^, and fluctuating interfaces either 

in the steady states or approaching steady states starting from random initial 
configurations 0. Typically one is interested in persistence, i.e., the probability 
Po{t) that at a fixed point in space, the quantity sign[(j){x,t) — {(j){x,t))] (where 
(j){x, t) is a fiuctuating field e.g., the spin field in the Ising model or the height of a 
fiuctuating interface) does not change upto time t. In all the examples mentioned 
above, this probability decays as a power law, Po(t) ~ t~^°, where the exponent 
60 is nontrivial. This nontriviality is due to the fact that the effective stochastic 
process in time at a fixed point in space becomes non-Markovian due to the cou- 
pling to the neighbours. For a non-Markovian process, calculation of any history 
dependent quantity such as persistence is extremely hard barring a few special 
cases [p^ , 14]. The exponent 6q has also been measured in a recent experiment on 



a liquid crystal system which has the same dynamics as the T = Ising model 
in 2-d|T0[. The experimental value was in good agreement with the analytical 
prediction of ^0 in 2-d Ising model[^. The exponent 6q has also been measured 
in a recent experiment on two dimensional soap froth [ITTI. 

In the above, one studied the persistence of a single spin (e.g., in the Ising 
or Potts model) of the initial random configuration. A natural generalization of 
this would be to study the persistence of a pattern, and not just a single spin, 
present in the initial configuration. Persistent patterns are quite abundant in 
nature. Examples include persistent eddies and vortices in turbulence, the great 
red spot of Jupiter or certain patterns of stock prices in financial markets. Another 
example is the so called activity-centered pattern in a self-organized system such 



as an interface in a random medium [Ol and also in certain models of evolution IIT? 



A natural question then is: what is the probability that a given pattern survives 
upto time t? 

Such persistent patterns exist also in phase ordering systems such as the q- 
state Potts model. For example, one such pattern is an original domain of a 
specific colour present in the random initial configuration of the Potts model. 
One can then ask: what is the survival probability of such a domain upto time t? 
This quantity for the 1-d Potts model has recently been studied by Krapivsky and 
Ben-Naim[Bn]. However this can be a more general question for any fiuctuating 
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field such as the solution of diffusion equation with random initial configuration 
or a fluctuating interface approaching the steady state. In such examples, a 
domain would be a connected set of points where the sign of the fluctuating field 
is positive (or negative). Another example of "pattern" persistence would be to 
study the probability that two adjacent domains in the initial configuration both 
survive upto time t. In this paper we develop a general framework to study the 
persistence of patterns of a fiuctuting field and discuss a few examples in detail 
where explicit results can be obtained. 

The general framework to study some of these "pattern" persistence problems 
consists of monitoring the motion of an external test particle launched in the 
fluctuating field. The dynamics of the test particle is suitably chosen so that the 
particle evaluates where the specific pattern of the fluctuating field is and moves 
there. The persistence of the pattern is then precisely the survival probability of 
the test particle. This led us naturally to study a more general "persistence of a 
mobile particle (P) in a field (F)" problem (henceforth the PF problem), special 
cases of which correspond to the "pattern" persistence in the underlying field. In 
this paper we study in detail a few examples of this general PF problem and find 
very rich, though often nonuniversal, behaviour. 

The general PF problem can be defined as follows. Let us consider a field 
(f){x, t) that fluctuates in both space and time. For example, the field 0(x, t) may 
be the solution of the simple diffusion equation 9^0 = V^0, or the order parameter 
profile of the Ising model undergoing T = coarsening dynamics or the height 
profile of a fluctuating interface or may even be spatially uncorrelated Brownian 
signal dt(j) = rj where ri{x,t) is spatially and temporally uncorrelated Gaussian 
white noise. A test particle is launched at an arbitrary initial point at time t = 0. 
The particle moves according to some prescribed deterministic or stochastic rules 
which in general depend on the local field profile. We now ask the question: given 
the dynamics of the particle, what is the probability P{t) that the field seen by 
the particle at its own location does not change sign upto time t? 

The survival probability of a mobile particle in a field has been studied before 
in the context of heterogeneous reaction-diffusion systems]^, ^ where the test 
particle was an external impurity diffusing through a homogeneous background. 
These studies dealt with a special case of our general PF problem, namely when 
the test particle is a simple random walker and the field is the coarsening "colour" 
field of the g-state Potts model at T = 0[|, |]. Amongst other studies on the 
similar line was the computation of the trapping time distribution of a diffusing 



tracer particle on a solid-on-solid surface|12|. In this paper, we extend these 
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studies to several other examples arising naturally in the context of our general 
PF problem. 

In the PF problem, one needs to specify the dynamics of the field as well 
as that of the test particle. For the field, we will consider three different cases, 
namely when the fields are: (i) the solution of diffusion equation, (ii) the spin 
profile of the Ising model or in general the "colour" profile of the g-state Potts 
model undergoing T = coarsening dynamics and (iii) spatially uncorrelatcd 
Brownian signals. The motion of the test particle, in general, consists of two 
separate moves. In the first part of the motion the particle sees the local field 
profile and then adopts a strategy to move in such a way so that it can live 
longer. This is the "adaptive" part of the motion that depends on the local 
field profile. In addition to this "adaptive" move the particle in general may be 
subjected to an external noise which consitutes the second part of the motion. 
This "noisy" part consists of Brownian moves of the particle that is independent of 
the local field profile. While the motion consisting of both "adaptive" and "noisy" 
moves is more general, for the sake of simplicity we will restrict ourselves to two 
separate cases when the motion is cither purely "adaptive" or purely "noisy" . In 
the first case, the survival probabihty is larger than the "static" case (when the 
particle is at rest) and in the second case, it is smaller. In cases where survival 
probabilities both in the mobile and the static cases decay as power laws, the 
exponent inequality 6ad ^ ^ (^rw holds where 6ad and 6rw are the persistence 
exponents associated with the "adaptive" and the "random walk" motion of the 
test particle. The relationship between the survival probability of the test particle 
and the persistence of "patterns" in the underlying field is established wherever 
possible. In all these studies, we will restrict ourselves to one dimension although 
in most cases the generahzations to higher dimensions is quite straightforward. 

The paper is organized as follows. In section- 1, we consider the "adaptive" 
motion of the test particle. Three cases of the fields are considered. By suitably 
choosing the adaptive strategy, the survival probability of the particle is related to 
that of a "pattern" (original domain in the Potts model or diffusion equation). An 
exact solution is presented for the special case when the motion of the "adaptive" 
particle is "directed" and the field is spatially uncorrelatcd Brownian signals. In 
section-2, we consider the "noisy" motion of the test particle where it performs 
simple "random walk" . Approximate analytical methods are developed to calcu- 
late the exponent characeterizing the power law decay of the survival probability 
of the particle. In section-3, a special case of the survival of the the diffusive 
test particle in the 1-d Potts model is shown to be related to the persistence of 
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an initial pattern, namely two adjacent original domains. This is also related to 
the fraction of uncollided domain walls at time t in the 1-d Potts model which 
we study both by numerical and analytical methods. Finally, we conclude with 
a summary, some suggestions for future directions and possible experimental re- 
alizations or tests of our results. Some details about a variational calculation are 
presented in the Appendix. 

II. Adaptive Motion of the Test Particle 

In this section we consider the adaptive motion of the test particle in which the 
particle adopts a strategy to move in a way such that it can live longer. The 
model and the strategy is as follows. Consider a lattice with periodic boundary 
conditions for convenience. The field (f){i,t) evolves with time according to some 
prescribed dynamics. A test particle is launched at t = 0, at an arbitrary site say 
the origin. Let us assume that at t = 0, the sign of the field (p at the origin is 
positive (or negative). As time changes, the field (f){i,t) changes. As long as the 
sign of the field at the origin is positive (or negative), the particle does not move. 
When the sign changes at the origin, the particle looks for a nearest neighbour 
where the field is positive (or negative). If it finds such a neighbour it goes there. 
In case there are more than one neighbours with positive (or negative) fields, it 
chooses one of them at random. Then it waits there till the sign of at the new 
site changes and then again it moves to one of its current neighbours and so on. 
If at some stage, the sign changes at the particle's current site and it fails to find 
a neighbour with positive (or negative) field, then it dies. Then we ask: What 
is the probability Pad{t) that the particle survives upto time t. Note that if the 
particle did not move at all and stayed put at one site only, then the survival 
probability Poit) is the usual "static persistence". 

We first consider the case when the field 0(x, t) is the solution of the simple 
diffusion equation dt4> = V^0, starting from a random initial configuration of 
0. For this simple field, analytical computation of even the static persistence 
Po{t) turned out to be quite nontrivial^. In 1-d, it was found that Po(t) ~ t^^° 
for large t where 6q = 0.1207 ± .00050]. The reason for the nontriviality once 
again can be traced to the fact that the effective Gaussian process that a static 
particle sees in time is non-Markovian. Nevertheless, an "independent interval 
approximation" (IIA) was developed in|Q, which produced analytical predictions 
for 9q for all dimensions that were extremely accurate. 

In case of the moving particle with the adaptive strategy, we carried out a 
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numerical simulation. The results are presented in Fig. [1]. We find, Pad{t) ~ 
t-^-d for large t, where Bad = 0.091 ±.002, compared to = 0.12 ±.001. Thus the 
strategy adapted by the particle is "successful" in the sense that the exponent, 
and not just the amplitude, characterizing the decay of persistence decreases. 
Thus in the renormalization group language, the adaptive strategy is a relevant 
perturbation. 

It is clear from above that the strategy that the particle adopts is basically 
to move towards the local maximum (or minimum) of the underlying field if 
the initial sign of the field that the particle sees is positive (or negative). By 
symmetry of the initial condition the survival probability of the particle does 
not depend on the initial sign of the field that it sees. Hence without any loss of 
generality it is sufficient to consider the case when the particle moves only towards 
local maximum. This observation may be used to develop a possible continuum 
approach to this problem. Let R{t) denote the position of the particle measured 
from a fixed point in space and x denote the co-ordinate of an arbitrary point in 
space measured from the location of the particle. Then, the effective field ilj{x,t) 
as seen by the particle is given by, ip{x,t) = (^{x ± R{t),t). Note that x = 
denotes the position of the particle. Given that (p satisfies diffusion equation, the 
equation of motion of ip{x^t) is given by, dtip = V'^ip ± R{t)dxip- We now model 
the adaptive strategy (namely that the particle tries to move towards the local 
maximum) by assuming that the velocity of the particle is proportional to the 
local slope of the field that the particle sees, i.e., R = Xip'{0,t), where ip' denotes 
the derivative with respect to x and A is a constant. The implication of this 
assumption is clear. If the local slope is positive, the particle moves to the right 
and if the local slope is negative, the particle moves to the left. Thus the particle 
always tries to move towards the local maximum. So, the equation satisfied by 
?/'(a;, t) is 

dtijj = V^ijj + Xijj'ix, t)^'{0, t) (2.1) 

which is a non-local and nonlinear KPZ type of equation. Then the adaptive 
persistence in this formulation is the probability that the local field V'(0,t) does 
not change sign upto time t. While a continuum formulation does not make it 
easier to compute the adaptive persistence exponent, but it relates the problem 
to the more familiar problem of persistence of fiuctuating interfaces^]. However 
we will not study this equation any further in the present paper and will defer its 
discussion for future [|19|. 

We now turn to the case when the field is the "colour" field of the g-state Potts 
model undergoing T = temperature dynamics. At each site of a lattice, the 
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field can take q possible colours. One starts from a random initial configuration 
of colours. A site is chosen at random and its colour is changed to one of its 
neighbours. This is how the colour field evolves. A test particle is launched as 
usual and it waits at its initial site till the colour of that site changes. Then it 
tries to find a neighbour with the same colour and if succeeds it goes to that 
neighbouring site. When it does not find any neighbour of its own colour, the 
particle dies. Then the question is as before, what is the probability Padif) that 
the particle survives upto time t. The corresponding "static persistence" exponent 
has been calculated exactly for all q recently^. 

Our job for calculating the adaptive persistence for the g-state Potts model 
is simplified by making the observation that the test particle survives as long as 
the original domain of the Potts model that contained the test particle at t = 
survives. Thus the adaptive persistence is precisely the survival probability of an 
"original domain" that has recently been studied both numerically and analyti- 
cally within an independent interval approximation (IIA) for all q by Krapivsky 
and Ben-NaimpO|. In fact, even for the diffusion equation, adaptive persistence 
is also the survival probability of an original domain. However, the IIA approx- 
imation developed in|2^ for the Potts model can not be easily extended to the 
diffusion equation for the following reason. The evolution of the Potts model is 
particularly simple in terms of the domain walls where the field changes colour in 
space. These domain walls perform independent random walks (the rates of which 
do not depend on the local spins) and when two walls meet, they either annihilate 
(with probability l/(g — 1)) or aggregate (with probability (g — l)/(g — 2)). There- 
fore, for the Potts model it is quite simple to write down an evolution equation for 
P{n, m, t) (probability that a domain of length n contains m original domains) 
within IIA[^ and thereby calculate the domain survival probability. However, 
writing down a similar evolution equation for the diffusion equation does not 
seem to be easy as the domain walls in the diffusion equation (i.e., the zeroes of 
the diffusive field) undergo complicated motion that depends upon the local field 
profile. 

However we want to stress that the concept of "adaptive persistence" is more 
general than just being equivalent to the persistence of a pattern, e.g., the domain 
survival probability in case of Potts model or diffusion equation. In fact, it does 
not necessarily require that the evolving field has domain structures coarsening in 
time. For example, even in the simplest case where at each site of the lattice there 
is an independent Brownian signal (completely uncorrelated spatially), one can 
define the "adaptive persistence" and has a nontrivial exponent as we will show 
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below via an exact solution. In cases when the evolving field has a structure of 
coarsening domains, the adaptive persistence is equivalent to the domain survival 
probability. 

We now turn to the last case for which we present an exact solution of the 
adaptive persistence. In this case, the field is spatially uncorrected Brownian sig- 
nal at each lattice site, dt(t){i,t) = r]{i,t) where ?7(i,t)'s are spatially uncorrected 
Gaussian white noises with zero mean and {v{ht)vUyt')) = ^ij^i't ~ t')- A test 
particle is launched as usual at t = at the origin. Let us assume that the signal 
at the origin at t = is positive. The particle does not move as long as the sign of 
the signal at the origin does not change sign. When it does, the particle either 
dies with probability (1 —p) or survives with probability p and then tries to jump 
to its neighbour on the right hand side. If the sign of the signal at that neighbour 
is positive at the time of jumping, the particle stays there till the signal is positive 
there and so on. Finally, if it does not find a right neighbour with a positive signal 
at the moment of jumping, the particle dies. Note the two new aspects in this 
problem from before. A survival factor p is introduced, p = 1 is the fully adapted 
case considered earlier for the diffusion equation or the Potts model, p = will 
correspond to the "static persistence". The second aspect is that the motion of 
the particle is directed as opposed to the undirected case considered for the diffu- 
sion equation or the Potts model. This assumption of directedness turns out to 
be important for exact solution and thus serves as an useful exactly solvable toy 
model of adaptive persistence. It turns out, as we show below, that the exponent 
6ad can be calculated exactly and depends on the parameter p continuously. 

Let Poit.to) denote the probability that a Brownian signal at a given site 
does not change sign from time to to time t. This is the usual static persistence, 
which can be computed exactly since it is a Markovian process and is given by, 
-Po(^0! ^) = f sin~^(^^^y=^)|T^. Let -Fo(to, t) denote the probability that the signal 
crosses zero for the first time at time t. Then clearly, Fo(to,t) = —dP^/dt. Then 
the probability of no zero crossing P{to,t) for the adaptive particle is given by 
the convolution, 

P{to, t) = Poito, t) + (p/2)Fo * Po + (p/2)'Fo * Fo * + • • • , (2.2) 

where Fo*Pq = J Fo(to,ti)dtiPo(ti,t) and so on. The first term is the probability 
that the sign of the signal at the starting site did not change upto time t and 
hence the particle did not jump at all. The second term denotes the probability 
that the particle jumped once. The parameter p is the survival factor and 1/2 
is the probability that the sign of the signal of the right neighbour (to which 
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the particle jumps) is positive. The third term denotes the probabihty that the 
particle jumped twice and so on. 

To perform the convoluted integrals in Eq. ( f^.2| ) , we make a change of vari- 
able, Tj = loglti/to). In this new variable, the effective process that the particle 
sees, though not Gaussian, becomes stationary. One can then use the Laplace 
transforms to solve Eq. ( |2.2D . Let -P(s), Po{s) and Fq{s) denote the respective 
Laplace transforms of -P(T), Po{T) and Fo{T). Then, by taking the Laplace 
transform of Eq. (|2^), one gets. 

Pis) = -^^k)—. (2.3) 

Using the realtion, Fq(T) = —dPo/dT, one further gets, Fq{s) = 1 — sPo{s). This 
enables us to write P{s) entirely in terms of Po{s), 



p(s) = . (2.4) 

We expect P{to, t) to decay as t"^*"* for large t. This means that in the variable 
T = logit/to), P(T) ~ exp(— 6'arfT) for large T. This implies that the Laplace 
transform P{s) will have a pole at s = —Bad, i-e., the denominator of the right 
hand side of Eq. (|2.4| ) will have a zero at s = —9ad- Using Po{T) = ^ sin~^(e~"^/^), 
one therefore sees that that the exponent 9ad is given by the positive root of 

l_P_P^ r sin\e-^'^)e'^'^^dT = 0. (2.5) 



2 vr Jo 

This integration can be evaluated by parts and one finally gets, 

B[l/2, 1/2 - 9ad] = — (2.6) 
P 

where B[m,n] is the usual Beta function. It is clear that in the limit, p 0, 
one recovers the usual "static" persistence exponent = 1/2. For the fully 
adapted model [p = 1), we get 6ad = 0.3005681 . . . which agrees very well with 
our numerical simulations. It is also clear that for any nonzero p, the adaptive 
exponent 9ad < 6'o = 1/2 as expected. 

The reason that the exponent 6ad is exactly soluble for the directed case is 
that the effective process seen by the test particle is Markovian. In the undirected 
case, this is not so because the particle can jump back to a site already visited 
before and therefore the probability that the signal is positive there at the time 
of current jumping is conditioned by the fact that the signal had crossed zero 
there at some earlier time. Therefore, it is difficult to compute 6ad exactly for the 
undirected case. 
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III. Noisy Motion of the Test Particle 



In the previous section we considered the adaptive motion of the test external 
noise. In this section we consider the other case when there is only noise and 
no adaption. In this case the particle just moves randomly through the medium 
in which a field (f){x, t) is evolving according to its own prescribed dynamics. As 
before, the test particle is launched at the origin at t = where the sign of the 
field is positive (say) at t = 0. The particle then performs a Brownian motion, 
R = ri{t), where //(t) is a Gaussian white noise with zero mean and the correlator, 
{ri(t)ri{t')) = 2Dp6(t — t'). Here R{t) denotes the position of the particle from 
some fixed reference point. The motion of the field and that of the particle are 
completely uncorrelated. The particle dies when the field that it sees changes 
sign. Then we ask: what is the probability that the particle survives upto time 
tl 

We first consider the case when the field t) is evolving according to diffu- 
sion equation, dtcf) = DfV'^(j). We performed numerical simulation to compute the 
survival probability of the test particle. This probability decays as a power law, 
P{t) ~ t^^'* for large t where the exponent 6d is found to depend continuously on 
the ratio of the two diffusion constants, c = Dp/ Df. The results are presented in 
the second column of Table (1) under the heading 6d{MC). It is clear from this 
table that for any nonzero c, 9d{c) > Oq where 6*0 is the corresponding "static" 
persistence exponent, i.e., when c = 0. This continuous non-universal dependence 
of 6d on c is, however, not very surprising for the following reason. The test par- 
ticle dies whenever it crosses any "zero" of the field (j){x,t). Thus two zeros of 
the field on either side of the test particle act like two boundary walls. But these 
walls are not static. They themselves are moving as the field t) is evolving 
in time. In fact, since the typical distance between zeroes of the diffusing field 
increases as ^/i, the walls bounding the test particle are therefore diffusing as 
This particular case is known to be marginal [p^ in the sense that the exponent, 
characterizing the power law decay of survival probibility of the particle, depends 
continuously on the ratio of the diffusion constants of the particle and the walls. 

While this explains qualitatively why 9d depends continuously on c, it does 
not give any quantitative estimate of the exponent. To make progress in that 
direction, we proceed as follows. Let x be the coordinate of an arbitrary point 
in space measured from the rest frame of the particle. Then the field il){x,t) = 
(f){x + R{t),t) as seen by the particle evolves as, 

|^...VV.f,W. (3.1) 
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The field ijj and the noise rj are completely uncorrelated. Note that for a given 
realization of the noise process {Ti{t)}, the process ip{x, t) at a fixed x as a function 
of t is a Gaussian process. However, when the distribution of 77 (t) is also taken 
into consideration, ip{x,t) at a fixed x no longer has a Gaussian distribution due 
to the multiplicative nature of the noise in Eq. (|3.1| ). It is nevertheless useful to 
calculate the the two time correlator, C{t',t) = (^/^(O, t')^/'(0, t)) that characterizes 
the temporal process at the location of the particle, i.e., at x = 0. This can be 
easily performed in the k space where the solution is given by, 

= V^(fc,0)e-^^'^V^/o^(*')'^*'. (3.2) 

We then compute C{t',t) = J dk{ilj{—k,t')ilj{k,t)) where the () is done over both 
the initial conditions of and the history of the noise rj. The initial consition is 
taken to random, so that {iplk, 0)ilj{—k, 0)) = A, where A is a constant. Since 
the noise i] is Gaussian white noise, we use the property, {exp(ik J^, r]{ti)dti) = 
exp[—Dpk'^\t — 1'\]. It is then easy to see that in d = 1, 

Cit\t) = V , ^ (3.3) 

jDf{t + t') + ap\t-t'\ 



The normalized autocorrelator, f{t\ t) = C{t', t) / yC{t', t')C{t, t) when expressed 
in terms of the variable T = log{t), becomes stationary, i.e., only a function of 
the time difference |T — T'|. Denoting, for convenience, this time difference by T, 
one finds that the stationary autocorrelator f(T) is given by, 

f{T) = I ^ (3.4) 

/cosh(T/2) + csinh(|T|/2) 



where c = Dp/Df is the ratio of the two diffusion constants. 

Note that for c = 0, f(T) reduces to the "static" autocorrelator /o(T) = 
[cosh(T/2)]"i/2||, However, there is an important difference between c = and 
c 7^ cases. For c = 0, the stochastic stationary process, whose correlator is 
given by foiT), is Gaussian. However, for c 7^ 0, while the process is still sta- 
tionary in the variable T, it is non-Gaussian as mentioned earlier. This is evident 
from Eq. (p.2|) , since ^jik^t) is a product of two random variables, ip{k,Q) and 
exp[iA; / rj{t')dt'] and hence is not clearly Gaussian, even though both 'ipik, 0) and 
Tjit) are individually Gaussian. Therefore, the approximate method developed^] 
for the case c = 0, to compute the asymptotic distribution of the interval between 
successive zero crossings, PoiT) ~ exp(— ^o^) ~ for large t, can not be simply 
extended to the c 7^ case. 
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It is nevertheless useful to calculate 6d by approximating the process by a 
Gaussian process having the same two-point correlator for two reasons. First, 
comparision between 6^ obtained numerically for the actual process and that 
obtained for the corresponding Gaussian process will tell us how important are 
the non Gaussian effects. Secondly, there have been some recent development s |1T5[] 
in approximate analytical calculations of the exponent for Gaussian stationary 
processes which one can use in the present context. Therefore, in the following, 
our strategy would be to to estimate the exponent 0(i{G) (that characterizes the 
exponential decay of the distribution of intervals between successive zero crossings, 
P(T) ~ (>-^d{G)T large T) for the Gaussian process with the correlator /(T) 
as in Eq. ( |3.4| ) and then compare it with the 6d for the actual process. 

We first present the numerical results for 9fi{G) for the Gaussian process with 
the correlator as in Eq. ( |3.4| ). This is done by constructing a time series having 
the same corelation function. It is most conveninent to work in the frequency 
domain rather than time domain. Details of this simulation procedure can be 
found in Ref.[^. The results of 6d{G) for different values of c are presented in the 
third column of Table (1). By comparing column (2) and (3), it is evident that the 
non-Gaussian effects are indeed quite small, and Gaussian approximation seems 
to be quite good. 

However, exact analytical calculation of 6d{G) even for a Gaussian stationary 
process with a general corelator /(T) is difficult and remains an unsolved problem 
for many years|TB|, Exact results are known only in a few special cases| 



One such case is when /(T) = e"^"!^! for all T. In this case, the Gaussian process 
is a Markov process and one can show exactly that P{T) ~ e~^^ for large T where 
e = Ao. By looking at /(T) in Eq. (gj), we see that for c = 1, /(T) = e"!^!/^ 
for all T. Therefore, 6d{G) = 1/4 for c = 1. For c close 1, say c = 1 + e, one 
can use a perturbation theory that has recently been developed to calculate 6d for 
processes that are close to a Markov process ITS], pT| . According to this theory, if 
f(T) = exp(— A|T|) + e/i(T), where e is small, then the exponent 9d, to order e, 
can be most easily expressed as, 

OdiG) = A[l - e— / /i(r)[l - exp(-2Ar)]-3/2rfr]. (3.5) 



In our case, expanding /(T) in Eq. (|3.4|) around c = 1, we get A = 1/4 and 
/i(T) = — i sinh(T/2) exp(— 3T/4). Using this in Eq. ( |3.5| ) and performing the 
integration, we get, to order e, 

^.(G') = i + 6A (3.6) 
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This perturbation theory may not give good estimates for 9d{G) when c is far 
away from 1. However, one can use a variational estimate for 6^(0) for general c. 
This variational method was recently developed in Ref. by mapping the zero 
crossing problem to that of the evaluation of ground state energy of a correspond- 
ing quantum problem. This method was used[|l5| to approximately calculate the 
"static" persistence exponent for the Ising model both in 1 and 2 dimensions. The 
results were in good agreement [|15|] with the exact result in l-d[0 and numerical 
simulations ll^] as well as direct experiment ||1CI|| in 2-d. This method works for 
class-1 Gaussian stationary processes, i.e., when /(T) = 1 — a\T\ + . . . for small 
T. Since in our present case, /(T) in Eq. (|3.4| ) is class-1 for any nonzero c, one 
can estimate 9d{G) by using this variational method. This method gives two esti- 
mates 6 max and Ovar for the exponent 6d{G). While 9max is a strict rigorous upper 
bound for 9d{G), Oyar gives the best variational estmate. The salient features of 
the variational method and the expressions for Omax and Oyar are given in the 
Appendix. 

For the process being considered with the correlator as in Eq. ( p.4|) , estimates 
Ovar and Omax are presented respectively in column (4) and (5) of Table (1) for 
different values of the parameter c. Comparing these with columns (2) and (3), it 
is clear that the variational approximation gets progressively worse as c increases. 
A visual summary of these different measures of the exponent is given in Fig. (2). 

To summarize, we find that the exponent 9d depends continuously on the ratio 
c of the diffusion constants. For any arbitrary nonzero c, 6d > Oq, where 6q is the 
corresponding static persistence exponent. Non-Gaussian effects are found to be 
quite small. 

We now turn to the case when the fluctuating field is the spin field of the Ising 
model or in general, the "colour" field of the g-state Potts model undergoing 
T = coarsening dynamics starting from a random initial configuration. The 
tracer particle, once again, moves diffusively through the medium with a diffusion 
constant Dp and whenever the field that the particle sees at its own location 
changes sign, the particle dies. As before, one is interested in calculating the 
survival probability P(t) of the tracer particle. 

In the g-state Potts model at T = 0, the domain walls perform random walk 
and whenever two domain walls meet, they either annihilate each other with 
probability l/(g — 1) or coagulate to form a single wall with probability {q — 
2)/(9 — 1)0) iH- The tracer particle dies whenever it crosses path with any 
domain wall. One expects that the survival probability of the tracer particle will 
decay as Pit) ~ t~^p for large t. The exponent 9p{q, c) is expected, as in the 
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diffusion case, to depend continuously on q and the ratio c = Dp/ Df where Dp is 
the diffusion constant of the test particle and -D/ is that of the domain walls. 

This problem has been studied in some detail before]^ ^. Let us just summa- 
rize here the main results that are already known. In the limit Q — > oo, 9p(^oo, c) 
can be computed exactly by noting that only two domain walls on either side 
of the tracer particle actually matter for the calculation of -P(t)[^. One finds 
exactly, 9p{oo, c) = 7r/{2cos~"^[c/(l + c)]}[|^, |^. Note that for c = 0, this reduces 
to the static persistence exponent, Op{oo, 0) = 10. In the Ising limit, q = 2, how- 
ever, there is no exact result for general c. Only exact result is available for c = 0, 
9p{2,0) = 3/80. For general c and g = 2, a mean field Smoluchowski type ap- 
proached was developed[|], whose predictions, 9p{2, c) = ^^(1-1- c)/8 were in good 
agreement with numerical simulations^. However, this Smoluchowski approach, 
when extended to large q limit, differed substantially from the exact q —* oo 
result. Finally a perturbation theory has been developed recently by Monthus^, 
and the exponent 9p{q, c) has been determined at first order perturbation in (<? — 1) 
for arbitrary c and at first order in c for arbitrary q. 

We will not study this exponent 9p{q,c) in its generality any further in this 
section. However, in the next section, we will study in some detail the special 
case c = 1 , as it turns out to be a particularly interesting case from the point of 
view of the "pattern" persistence problem of the Potts model. 

IV. Persistence of a Specific Pattern in the 1-d 
Potts Model 

In section-2, we showed that the survival probability of an "adaptive" test particle 
is related to the peristence of a specific pattern namely an original domain in the 
T = dynamics of the Potts model. In this section, we show that the survival 
probability of a "noisy" or "diffusive" test particle (studied in section-3) is also 
related to the persistence of yet another pattern in the 1-d Potts model namely 
the survival upto time t of two adjacent original domains present in the initial 
configuration. 

Let us consider the zero temperature coarsening dynamics of the g-state Potts 
model starting from a random initial configuration. In an infinitesimal time in- 
terval dt, each spin changes its colour to that of one of its neighbours selected at 
random. This dynamics can be equivalently formulated in terms of the motions of 
domain walls. The domain walls perform independent random walk and whenever 



14 



two walls meet, they either annihilate with probability — 1) or aggregate to 
become a single wall with probability (g — 2)/(g — 1)0, |2^. The "static" per- 
sistence then is the probability that a fixed point in space is not traversed by 
any domain wall. However a somewhat more natural quantity, in this domain 
wall representation, is the probability Pi{t) that a given domain wall remains 
uncollided upto time t. A little thought shows that this is precisely the survival 
probability of two adjacent domains present in the initial configuration. We show 
below that -Pi(t) ~ t~^^ for large t where Oi{q) is a g-dependent exponent that is 
not obviously related to any other known exponent via scaling relations. 

In section-2, we considered the exponent 6'p(g, c) characterizing the decay of the 
survival probability of a diffusing test particle in the background of the diffusing 
domain walls of the Potts model. The parameter c is the ratio of the diffusion 
constant of the test particle to that of the domain walls. Let us consider the case 
c = 1. In this case, the test particle can not be distinguished from the other 
diffusing domain walls. Since the test particle dies whenever any other domain 
wall touches it, it is clear that the survival probability of the test particle for c = 1 
is precisely the fraction of uncollided domain walls in the Potts model and hence 

Clearly 9i can be exactly determined in the two limits, q = 2 and g — * oo. For 
q = 2, since there is only annihilation upon contact between domain walls, the 
fraction of uncollided walls is the same as the density of domain walls that decays 
as ~ t~^/^ for large t, and hence 6*1(2) = 1/2. In the g — 00 limit, by putting c = 1 
in the exact formula, 9p{oo,c) = 7r/{2cos^^[c/(l + c)]}|^, |^, one gets 9i = 3/2. 
For intermediate values of g, we present numerical results in the column (2) of 
Table (2). It is clear that the exponent 6*1 (g) increases monotonically with g. 

The exponent 9i can be quite easily computed within mean field theory. Let 
N{t) and Ni{t) denote respectively the total density of domain walls and density 
of "uncollided" walls at time t. Let Qi{t) denote the density of domains of size 
1 (here 1 is the lattice spacing and hence the smallest interval size). Then N{t) 
and Qi{t) are related via the exact relation, 

dt g — 1 

However there is no such simple exact relationship between Ni{t) and Qi{t). 
However, if one neglects correlations, it is easy to write such a relationship within 
mean field theory, 

^ = -2(^m-2^(l-^)Q, (4.2) 
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where the two terms on the right hand side are self explanatory. Eliminating 
Qi from the two equations above, we get, A''i ~ Using the result, 

N{t) ~ t^^/^, we finally obtain, 9i{q) = {q — l)/q, within mean field theory. 
While the meanfield answer is exact for g = 2, it gets worse as q increases as 
evident by comparision with Table (2). Presumably the mean field value forms a 
lower bound to the true exponent value, though we have not been able to prove 
it. 

However, one can obtain rigorous upper bounds to 6i{q) as follows. This 
can be done by generalizing the arguments used by DerridaQ to obtain upper 
bounds to the "static" persistence exponent 9o{q). The argument goes as follows. 
It was noted by MonthusP] that the problem of a diffusing tracer particle mov- 
ing amongst the domain walls of the Potts model can be mapped to a reaction 
diffusion problem where particles are generated from a source, diffuse around and 
aggregate upon contact. The only difference from the "static" case[^ was that 
the source is now "moving". In fact the source diffuses with the same diffusion 
constant as the tracer particle. It is then possible to write the survival probability 
Pi{t) asi, 

oo 

Pi(t) = 5:P(m,t)gi— (4.3) 
1 

where P{m, t) is the probability of having m particles in the corresponding re- 
action diffusion problem. Writing the above equation for q = q2 in the following 

way, a > 1], 

P (n ) °° 

= J2 P{m,t)q^"' = ;^P(m,t)[gr'"]'°s92/iog9i^ (4 4) 

92 rn=l 

and then using Jensen's inequality [(x") > (a;)" for a > 1 and x a positive random 
variable] as was used in the "static" case0, we immediately obtain the following 
inequality, 

^ife) < e,{q,)^^ (4.5) 
loggi 

for q2 > qi- For example, using the exact result, 6'i = 1/2 for g = 2 and the above 
inequality we get, 

^i(g) <logg/21og2. (4.6) 

For g = 3, this gives ^i(3) < 0.792481..., which should be compared with its 
numerical value 0.72 ± 0.005. For higher values of q, one can have a numerical 
estimate of a tighter upper bound of 0i{q) by using the numerical value of ^i(g — 1) 
in the inequality ( [4.5|) . This would be an improvement over the exact bound (f4.6|) 
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obtained by comparing with q = 2. These improved numerical bounds 9i{max) 
are presented in the column (3) of Table (2). 

It was pointed out by Monthus[^] that carrying out the same formalism that 
led to the exact determination of the static persistence exponent ^^o(q')[0] is not 
so straighforward to compute the exact value of 9p{q, c) for general c. However, 
one may hope that some special simplifications might occur for c = 1 leading to 
the exact computation of the exponent 9i (g) though we have not succeded yet in 
that direction. 



V. Summary and Conclusion 

In this paper, we have studied the persistence of some patterns present in the 
initial configuration of a fluctuating field. It was shown that some of these pattern 
persistence problems are related to the survival probability of a mobile particle 
launched into the field. By suitably adjusting the rules of the dynamics of the 
particle one can study persistence of different patterns in the underlying field. 
This led us to study the PF problem in general. Several special cases were studied 
in detail and new results were derived. 

It is clear from ours as well as other studies that there is a whole hierarchy 
of exponents associated with the decay of persistence of different patterns in 
the phase ordering systems. It is not clear at present whether these exponents 
are independent of each other or not. While these exponents do not depend on 
the details of the initial configuration (as long as it is short ranged), it is not 
clear whether they can be considered "universal" and if so, in what sense. For 
example, if the diffusion constant of any single domain wall of the Potts model 
changes slightly, then the exponent 9i characterizing the decay of survival of the 
wall (probability that it remains uncollided) also changes. Clearly in this respect 
the exponent is non- universal. So the important question that remains to be 
answered is: what are the criteria one should use to decide whether an exponent 
in the phase ordering dynamics is universal or not? 

One of the interesting extensions of the present work would be to study the PF 
problem when the fluctuating field is the height of an interface in or approaching 
the steady state. The "static" persistence for interfaces have been recently studied 
in some detaiip. Also anomalous diffusive behaviour of a tracer particle on a 



solid-on-solid surface was noted [1^ and was attributed to the temporary trapping 
or burial of the particle in the bulk of the crystal. Besides, given that sophisticated 
techniques using scannning tunneling microscope already exist for determining 
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temporal step fluctuations on crystal surfaces p5[], it is not unreasonable to hope 
that such techniques may be refined in future to measure the survival probabilities 
of the "static" as well as the "mobile" particle in a fluctuating interface. 

Finally it has recently been noted that "static" persistence exponent for the 
diffusion equation may possibly be measured in dense spin-polarized noble gases 
(3He and 129Xe) using NMR spectroscopy and imagingj^. In these systems the 
polarization acts like a diffusing field. With a slight modification these systems 
may possibly be used to measure the persistence of some patterns of the diffusive 
field as discussed in the present paper. 
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Appendix: The Variational Method 

It was shown in that the exponent 6^(0) is exactly the ground state energy 
difference, 6d{G) = Ei—Eq between two quantum problems, one with a hard wall 
at the origin and the other without the wall. The energy Eq (without wall) can 
be exactly determined. 



ir,og(^M. (5.i; 



where G{uj) = l/f{uj) and f{uj) is the Fourier transform of the correlation func- 
tion /(T) of the Gaussian stationary process normalized such that, f{uj) ~ cj"^ for 
large uj. The energy Ei (with a wall at the origin) is estimated variationally as it 
is hard to obtain exactly. For Class- 1 processes, one can use a harmonic oscillator 
with a wall as the trial state with the frequency of the oscillator as the varia- 
tional parameter. The variational energy Ei{ujq) is given by the expression [|15||, 



= c.o[^ + -(^ - 1) + - r dxi^^ -x'- mx)], (5.2) 

2 TC CUn vr Jo UJn 



where S{x) = J2'i^=i ncn/{x^ + 4n^) with c„'s given by, 

4 {2n)\ 
^" ~ vr22«(2n + l)!^n!(2n- 1) 



(5.3) 
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One then minimizes Ei{u!o) with respect to ujq, and uses this minimizing frequency 
ujmin to obtain the variational estimate, 9max = Ei{LJmin)—Eo{u!min) which also is a 
rigorous upper bound to the true exponent 9d{G) that characterizes the Gaussian 



process. However, as argued in |T5[, one can obtain a better estimate of 9diG) by 



(2) (2) 

using, 9var = Ei{uJmin) - -Eq (u^mm), wherc is given by 



^i '--™»l2 + 2^/„ - 1)1. (5.4) 



19 



References 



[1] B. Derrida, A.J. Bray and C. Godrechc, J. Phys. A 27, L357 (1994); D. 
Stauffer, J. Phys. A 27, 5029 (1994); B. Derrida, P.M.C. de Oliveira and D. 
Stauffer, Physica 224 A, 604 (1996); E. Ben-Nairn, L. Frachebourg and P.L. 
Krapivsky, Phys. Rev. E 53, 3078 (1996). 

[2] B. Derrida, V. Hakim and V. Pasquier, Phys. Rev. Lett. 75, 751 (1995), and 
J. Stat. Phys. 85, 763 (1996). 

[3] B. Derrida, J. Phys. A 28, 1481 (1995). 

[4] S.N. Majumdar, C. Sire, A.J. Bray and S.J. Cornell, Phys. Rev. Lett. 77, 
2867 (1996); B. Derrida, V. Hakim and R. Zcitak, ibid. 2871. 

[5] A. Watson, Science, 274, 919 (1996). 

[6] P.L. Krapivsky, E. Ben-Nairn and S. Redner, Phys. Rev. E 50, 2474 (1994). 

[7] J. Cardy, J. Phys. A 28, L19 (1995); E. Ben-Naim, Phys. Rev. E 53, 1566 

(1996) ; M. Howard, J. Phys. A 29, 3437 (1996); S.J. CorneU and A.J. Bray, 
Phys. Rev. E 54, 1153 (1996). 

[8] C. Monthus, Phys. Rev. E 54, 4844 (1996). 

[9] J. Krug, H. Kallabis, S.N. Majumdar, S.J. Cornell, A.J. Bray and C. Sire, 
to appear in Phys. Rev. E (1997), cond-mat 9704238. 

[10] B. Yurke, A.N. Pargelhs, S.N. Majumdar and C. Sire, Phys. Rev. E 56, R40 

(1997) . 

[11] W.Y. Tam, R. Zeitak, K.Y. Szeto and J. Stavans, Phys. Rev. Lett. 78, 1588 
(1997). 

[12] J. Krug and H.T. Dobbs, Phys. Rev. Lett. 76, 4096 (1996); P.C. Searson, R. 
Li and K. Sieradzki, Phys. Rev. Lett. 74, 1395 (1995). 

[13] D. Slepian, Bell Syst. Tech. J. 41, 463 (1962). 

[14] M. Kac, SIAM Review 4, 1 (1962); I.E. Blake and W.C. Lindsey, IEEE Trans. 
Information Theory 19, 295 (1973). 

[15] S.N. Majumdar and C. Sire, Phys. Rev. Lett. 77, 1420 (1996); C. Sire, S.N. 
Majumdar and A. Rudinger, unpublished. 



20 



[16] S. Krishnamurthy and M. Barma, Pliys. Rev. Lett. 76, 423 (1996). 

[17] P. Bak and K. Sneppen, Phys. Rev. Lett. 71, 4083 (1993). 

[18] S.J. Cornell, unpublished. 

[19] S.N. Majumdar and D. Dhar, unpublished. 

[20] P.L. Krapivsky and E. Ben-Nairn, cond-mat 9705155. 

[21] K. Oerding, S.J. Cornell, and A.J. Bray, Phys. Rev. E 56, R25 (1997). 

[22] P.L. Krapivsky and S. Redner, cond-mat 9603176; cond-mat 9502041, 
9505062. 

[23] S.N. Majumdar and D.A. Huse, Phys. Rev. E 52, 270 (1995); C. Sire and 
S.N. Majumdar, Phys. Rev. Lett. 74, 4321 (1995); Phys. Rev. E 52, 244 
(1995). 

[24] M.E. Fisher and M.P. Gelfand, J. Stat. Phys. 53, 175 (1988); D. bcn- 
Avraham, J. Chem. Phys. 88, 941 (1988); D. Considine and S. Redner, J. 
Phys. A 22, 1621 (1989). 

[25] N.C. Bartelt, T.L. Einstein and E.D. Williams, Surf. Sci. 312, 411 (1994). 

[26] Ronald L. Walsworth, private communication. 

[27] C.H. Tseng, S. Peled, L. Nascimben, E. Oteiza, R.L. Walsworth, and F.A. 
Jolesz, Harvard-Smithsonian preprint no: 4502, to appear in J. of Mag. Res., 
series B, 1997. 



21 



Tables 



c 


OdiMC) 


OdiG) 


Ovar 


^ TTICLX 


0.5 


0.20 ±0.01 


0.190 ±0.005 


0.189 


0.210 


1.0 


0.26 ±0.01 


0.25 = 1/4 


1/4 

/ 


1/4 


2.0 


0.35 ±0.01 


0.325 ± 0.005 


0.319 


0.350 


3.0 


0.42 ± 0.01 


0.389 ± 0.005 


0.363 


0.442 


4.0 


0.48 ±0.01 


0.439 ± 0.005 


0.396 


0.528 


5.0 


0.53 ±0.01 


0.496 ± 0.005 


0.422 


0.611 


6.0 


0.58 ±0.01 


0.527 ±0.005 


0.444 


0.688 


7.0 


0.62 ±0.01 


0.579 ±0.005 


0.463 


0.776 


8.0 


0.65 ±0.01 


0.628 ±0.005 


0.479 


0.839 


9.0 


0.69 ±0.01 


0.694 ±0.005 


0.493 


0.912 


10.0 


0.73 ±0.01 


0.723 ±0.005 


0.505 


0.984 



Table I. Estimates of the persistence exponents of a diffusing tracer particle 
through an external field evolving via diffusion equation. For different values of 
the ratio c = Dp/Df as shown in column (1), exponents are obtained (2) by 
direct Monte Carlo simulation of the process {9mc)) (3) by simulating a Gaussian 
stationary process with the correlator /(T) = [cosh(T/2)±csinh(|T|/2)]~-'^/^ {das) 
and (4) by using variational estimate Ovar for the above Gaussian stionary process 
with correlator /(T) (5) the rigourous upper bound O^ax for the above Gaussian 
process. Monte Carlo simulations were carried out on a periodic lattice of 100000 
sites and the results were averaged over 20 samples. 



22 



9 


9i 


6i{max) 


2 


1/2 


- 


3 


0.72 ± 0.01 


0.792481.. 


4 


0.86 ±0.01 


0.91 ±0.01 


5 


0.95 ±0.01 


1.00 ±0.01 


6 


1.04 ±0.01 


1.06 ±0.01 


50 


1.47 ±0.01 




oo 


3/2 





Table II. The exponent 9i that characterizes the asymptotic decay Pi{t) ~ t~^'^, 
probabihty that a domain wall remains uncoUided upto time t in the zero tem- 
perature dynamics of the g-state Potts model. Exact values of 9i are quoted for 
q = 2 and q ^ oo. For other values of q, 9i (column (2)) is estimated from Monte 
Carlo simulations on a periodic lattice of 75000 sites and results averaged over 20 
different initial conditions. The estimated upper bounds for 9i{q) (as explained 
in the text) are presented in column (3). 
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Figure captions 

Figure 1. Log-log plot of Monte Carlo simulations of the "adaptive" persistence, 
Pad{t) (plus symbols) and "static" persistence, Po{t) (cross) versus time t. The 
simulations were carried out on a periodic lattice of 100000 sites and results were 
averaged over 20 samples. The best fit to the straight lines gives the exponent 
values, Sad = 0.091 ± 0.002 and 60 = 0.12 ± 0.001 

Figure 2. A visual summary of the different measures of the exponent 9 a as 
given in Table 1. 
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